
capture log close
log using $log/EDS_10_robustness_inference, text replace

clear
clear matrix
clear mata
set mem 2g
set matsize 800
set more off
set maxvar 20000
set trace off
set varabbrev off


#d;
	  global cov "  _Isu_schule_1 _Isu_schule_2 _Isu_schule_3 
	_Isu_ausbila1 _Isu_ausbila2
	_Isu_kat_al_2 _Isu_kat_al_3 _Isu_kat_al_4 
	_Ichild_1 _Ichild_2
	_Isu_eintri_1 _Isu_eintri_2 _Isu_eintri_3 _Isu_eintri_4 _Isu_eintri_5
	_Isu_eintri_6 _Isu_eintri_7 _Isu_eintri_8 _Isu_eintri_9 _Isu_eintri_10 _Isu_eintri_11
	_Ibefstat_1 _Ibefstat_2 _Ibefstat_3 _Ibefstat_4
	_Idauloint__8 _Idauloint__9 _Idauloint__10 _Idauloint__11 _Idauloint__12 
	_Idauloint__13 _Idauloint__14
	female east_germany
    german migration_background 
    married  
	benefit_status benefit_level log_benefit_level
    lifetime_ue_age lifetime_empl_age last_wage last_wage_mis 
    region_ba_aloq_zivall gdp_real_pc region_ba_vacancies 
	sector_manufacturing sector_service  log_mobirate_avg5
	openness_st conscient_st extraversion_st neuroticism_st loc_st
	write_ger speak_ger write_eng speak_eng
	homeowner  
	   ";		

	   

		global outc"
		    search_dist50 
	       log_distance_km 
		    applications_total
		   applications_distant
		   applications_local
		   channel_active
		   ";
	  #d cr




		use $data/data_pair, clear

		joinby region_kkz using $data/koordinaten_kreise,unmatched(master)
		tab _merge
		drop _merge
		
		encode region_pair, gen(region_pair2)
		cap drop weight
		gen weight=1/fweight_reg

	
		
		foreach type in  $outc {
	 	
			
		
		qui reg `type'  log_intensity_alt  $cov  county_pair_*  [pw=1/fweight_reg], cluster(region)
		matrix b=e(b)'
		matrix n=e(N)
		matrix v==e(V)
		matrix se0=vecdiag(cholesky(diag(vecdiag(v[1..1,1..1]))))'
		
		boottest log_intensity_alt, reps(9999) cluster(region) nograph boot(wild)
	    matrix p = r(p)	  
		
		
		qui acreg `type' log_intensity_alt $cov [pw=weight], pfe1(region_pair2) spatial latitude(Breite) longitude(Länge) dist(25)
	    matrix b=e(b)'
		matrix b1 = b[1..1,1..1]
		matrix v==e(V)
		matrix se1=vecdiag(cholesky(diag(vecdiag(v[1..1,1..1]))))'
		
		qui acreg `type' log_intensity_alt $cov [pw=weight], pfe1(region_pair2) spatial latitude(Breite) longitude(Länge)  dist(50)
	    matrix b=e(b)'
		matrix v==e(V)
		matrix se2=vecdiag(cholesky(diag(vecdiag(v[1..1,1..1]))))'
		
		qui acreg `type' log_intensity_alt $cov [pw=weight], pfe1(region_pair2) spatial latitude(Breite) longitude(Länge)  dist(150)
	    matrix b=e(b)'
		matrix v==e(V)
		matrix se3=vecdiag(cholesky(diag(vecdiag(v[1..1,1..1]))))'
		
		matrix d_`type'=b,se0,se1,se2,se3,p
		
		
		
		
	
		}/*type*/
		
		
		
			#delimit;
		matrix d=
		d_search_dist50 \
		d_log_distance_km \
		d_applications_total \
		d_applications_distant \
		d_applications_local \
		d_channel_active
		;
	
		matrix rownames d1 = 
		"Applied to distant" 
		"Search radius" 
		"#applications total" 
		"#applications distant" 
		"#applications local" 
		"Active search channel"
		;
		#delimit cr
	matrix colname d = "Coef" "SE clustered" "SE spatial 25km" "SE spatial 50km"  "SE spatial 150km" "P-value WCB"
	estout matrix(d, fmt(4)) using $data/Table_A9_survey, style(tex) replace
	
	
	
log close	

